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We distinguish difTerent mechanisms for population inversion in flux qubits driven by dc+ac 
magnetic fields. We show that for driving amplitudes such that there are Landau-Zener-Stiickelberg 
intereferences, it is possible to have population inversion solely mediated by the environmental bath 
at long driving times. We study the effect of the resonant frequency Qp of the measuring circuit, 
flnding different regimes for the asymptotic population of the state of the flux qubit. By tailoring Q,p 
the degree of population inversion can be controlled. Our studies are based on realistic simulations 
of the device for the Josephson flux qubit using the Floquet-Born-Markov formalism. 
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Population inversion, where the highest populated 
state is an excited state, is among the most interesting 
phenomena in maser and laser physics •1;]. The usual 
mechanism to get population inversion requires driven 
quantum systems with three or more energy levels, since 
indirect transitions through a third energy level are in- 
volved [iHij. Here, we find that the device for the Joseph- 
son flux qubit (DJFQ) can have population inversion 
even in the case in which only two levels are involved 
in the driven dynamics. Furthermore when tailoring the 
structure of the environmental bath different mechanisms 
for population inversion can arise. 



Superconducting circuits with mesoscopic Josephson 
junctions are used as quantum bits and can behave as 
artificial atoms @. A well studied circuit is the DJFQ 
0, Q which, for millikelvin temperatures, exhibits quan- 
tized energies levels that are sensitive to an exernal mag- 
netic field. When driven by a dc-|-ac magnetic fiux, tran- 
sitions between energy levels at avoided crossings occur, 
and a rich structure of Landau-Zener-Stiickelberg (LZS) 
interferences |9| combined with multi-photon resonances 
is observed |10l4l2| . For large ac amplitudes diamond- 
like interference patterns are displayed, from which the 
energy level spectrum has been reconstructed 



LZS 



interferences have also been observed in charge qubits 
[ist . Rydberg atoms 



14| . ultracold molecular gases [1 



optical lattices [16| and single electron spins [17 1 



The DJFQ consists on a superconducting ring with 
three Josephson junctions Q enclosing a magnetic fiux 

= f^o i^o = h/2e) with phase differences (fi, (p2 
and ip3 — ~Lpi + Lp2 — 27r/. Two of the junctions have 
coupling energy, Ej, and capacitance, C, and the other 



has i?j_3 = aEj and C3 = aC. The hamiltonian is: 
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with pt = (<^i + <y32)/2, Pi = {pi - <y52)/2, Ec = eV2C 
and V{pt, 'Pi) — 2 + a — 2 cos ipt cos ipi — a cos{2tt f + 2pi) . 
The DJFQ has several discrete levels with eigenenergies 
Ei and eigenstates which depend on /, a and 77 = 
y/SEc/Ej. Typical experiments have a 0.6 — 0.9 and 
77 - 0.1-0.6 [iEllil. For a > 1/2 and |/ - 1/2| < 1, 
the potential V{ipt,pi) has the shape of a double- well 
with two minima along the pi direction. In this regime 
the system can be operated as a quantum bit B, @1 and 
approximated by a two- level system (TLS) 0,[l3l- When 
driven by a magnetic flux f{t) = fdc + /ac sin(a;i), the 
hamiltonian is time periodic 'H(t) = %{t + t), with r = 
27r/a;. In the Floquet formalism, that allows to treat 
periodic forces of arbitrary strength and frequency 
|20|, the solutions of the Schrodinger equation are of the 
form Yipit)) = e*^'3*/'*|$^(t)), where the Floquet states 
satisfy |ci>^(i)) = |$Mt+r)) = Y.k |1>^)e-'""*, and 
are eigenstates of the equation \H{t) — ihd / dt]\^ j3{t)) = 
e^|<i>^(i)), with ep the associated quasi-energy. 

Experimentally, the system is affected by the electro- 
magnetic enviroment that introduces decoherence and 
relaxation processes. A standard theoretical model to 
study enviromental effects is to cotiple the system with 



24- Following [24 



a bath of harmonic oscillators [21- 
we assume that the DJFQ is linearly coupled through pi 
to a bath that has spectral density J{ui) and is equi- 
librated at temperature T. For weak coupling (Born 
approximation) and fast bath relaxation (Markov ap- 
proximation), a Floquet-Born-Markov master equation 
for the reduced density matrix p in the Floquet basis, 
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Pap{t) = (<I>a(t)|p(t)|<J>/3(i)), can be obtained |2]| : 



dt 



£l3)Pal3{t) + La'p'aP Pa'l3'{t) 



a' 13' 



the time scale for dissipation is 



(2) 
larger 



the transition rates L^'fi'apit) can 



When 
than T 

be averaged over one period t '21], obtaining 

La'P'a(3 — J2kii^aa',k + Nj3i3'^k)^aa',k^l3'P,-k — 
J2,j {Spi3'N,^a',kTa7j,-k^j)a',k + Saa' Njjpi ^f; r],-k)] , 

with r^0^g = and N^^^g = iV(£„ - 

Sf) + qhuj), where N{s) = J(e/;i) (coth(e/2A:BT) - 1). 
This formalism has been extensively employed to study 
relaxation and decoherence for time dependent periodic 
evolutions in double- well potentials and in TLS l2l| - |23| . 
Here we use it to model the ac driven DJFQ, beyond 
the TLS approach, considering the full hamiltonian of 
Eq. JT]). We calculate the coefficients La'p'ap using 
the Floquet states |$^) and quasienergies £/3, obtained 
as in [20|, and then we integrate numerically Eq. 
obtaining Pajsit) as a function of t. 




FIG. 1: Population Pr as a function of the driving amplitude 
fac for t — lOOOr (continuous line). The dashed line shows the 
asymptotic (t — > oo) average population Pr. The inset shows 
the dependence of Pr on T for a driving amplitude fac ~ 
0.00245. The calculations were performed for fdc ~ 0.50151, 
to = 2n/T = 0.003-Ej/ft and an ohmic bath with 7 = 0.001 
at r = 20mK (for Ej/h f» 300GHz). Vertical lines separate 
regimes A, B, C, D described in the text. Black circle: fac 
corresponding to the inset, Fig.2(a) and Fig. 3; black square: 
fac corresponding to Fig. 2(b). 

We start by considering a DJFQ coupled to a bath with 
an ohmic spectral density J{uj) — "fUJ, which mimics an 
unstructurated electromagnetic environment that in the 
classical limit leads to white noise. We take 7 = 0.001, 
corresponding to weak dissipation as in [llj , and the bath 
at T = 0.0014£;j/A:i3 (- 20mK for Ej/h - 300GHz). 
The DJFQ parameters are a = 0.8 and 77 ~ 0.25. The 
static field is taken as fdc = 0.50151, and for the driving 
microwave field we choose uj — 0.003Ej /h ~ 900MHz 



and different amplitudes fac- The initial condition is the 
ground state {"^/q) of the static hamiltonian Hq = T-Lifdc)- 

Experimentally, the probability of having a state of 
positive or negative persistent current in the flux qubit 
is measured [Sj. The probability of a positive current 
measurement ( "right" side of the double- well potential) 
can be calculated as PR{t) — TxiJ\Rp{t)), with 11^ the 
operator that projects wave functions on the > sub- 
space For a static field fdc ^ 1/2, the ground state 
has Pfi K, 1. In Fig.l we show vs. faa calculated 
from the numerical solution of Eq. ^ after a driving time 
t = lOOOr, which is similar to the time scale of the ex- 
periments fllj, finding different regimes. (A) For small 
fac, Pr ^ 1, since the system is slightly perturbed from 
the ground state. (B) When further increasing fac, new 
regimes are found whenever the extreme driving ampli- 
tudes fdc ± fac reach an avoided crossing in the energy 
level spectrum {Ei{f))}. Since the slopes dEi/df are 
proportional to the average loop current of the DJFQ, 
there is an interference of "positive" and "negative" loop 
current states ("right" and "left" sides of the double- 
well potential) at avoided crossings, which results in LZS 
oscillations [9| in the dependence of Pr with fac- The 
first case is found in Fig.l when the avoided crossing at 
/ = 1/2 between the ground state level Eq and the first 
excited level Ei is reached, and the transfer of population 
to the El level lowers Pr. The DJFQ of [13, Ell has short 
decoherence times (i^ ^ 20ns > r ^ 1 — 10ns) and large 
relaxation times (t^ ^ 100/is). Due to this time scale 
separation {t < < t <s:_ tr), a. TLS model with classical 
noise, valid for t < tr, can explain the experimentally ob- 
served LZS oscillations of Pr vs. fac, where the minima 
correlate with the zeros of Bessel functions Jn{afac/oj) ( 
a a normalization constant) [l3|- The results of Fig.l 
for t = lOOOr are in agreement with this finding. (C) At 
higher values of fac more than two levels have to be con- 
sidered. When the avoided level crossing between Ei and 
E2 is reached by the driving amplitude the ground level 
Eq is repopulated due to fast E2 Eq 'intrawell' transi- 
tions (the corresponding states have the same sign of the 
average loop current) and Pr increases to values near 
1 ('cooling' effect, see [ll|). (D) For fac such that the 
symmetrically located (with respect to / = 1/2) avoided 
crossing between Ei and E2 is reached, there are new 
LZS oscillations. Furthermore the levels Eq and E^ are 
also involved in the dynamics (since the avoided cross- 
ings between E2 O E3 and Eq ^ Ei are traversed by 
the driving) and a more complex dependence of Pr with 
fac emerges. In this case, we find population inversion 
{Pr < 1/2), which is also observed experimentally in 1 1111 . 
A multilevel extension of the semiclassical model of [10[ 
can describe this behavior as well (25| . 

In Fig. 2 we show the explicit time dependence of Pa/sit) 
for two different values of fac- We find that the off- 
diagonal elements pa/s , go to zero exponentially for large 
t [see Fig. 2(a) and (b), left axis], with the decay time 
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FIG. 2: (color online) Off-diagonal elements of paja (left axis) 
and level population Pi (right axis) as a function of time (log) 
for fac = 0.50151, u = 0.003Ej/h, 7 = 0.001 and two differ- 
ent values of the driving amplitude, (a) fac = 0.00245 and 
(b) fac = 0.0191. 



corresponding to the decoherence time of the periodi- 
cally driven case. This shows that the density matrix is 
diagonal in the Floquet basis after full decoherence (26| , 
and the relaxation process is determined by the evolution 
of the diagonal Paa{i) for t > t^. Thus the asymptotic 
regime, t ^ tr, can be obtained from the non-trivial solu- 
tion of Lp/3,aaPi3i3 = after impossing dpaa/dt = in 
Eq.(2). In this way, we calculate Pn = limt^oo(^'fl'(i))T, 
averaged over one period r, which is shown in Fig.l 
(dashed line) as a function of fac- We see that, while 
for low and large values of /ac, it follows the finite time 
Pnit), for intermediate fac the asymptotic Pn is lower 
and can show population inversion. In order to under- 
stand the emergence of population inversion in the limit 
t tr, we calculate the time evolution of the popula- 
tion of the eigenstates of Hq, P^ = (*i|p(i)|*i). The 
results are shown in Fig. 2(a) and (b), right axis, where 
we distinguish two mechanisms of population inversion: 
(i) Bath-mediated population inversion: In Fig. 2(a) we 
show results for a value of fac such that only LZS in- 
terference between the two lowest energy levels occurs. 
There is a fast decay of pa/j (in the time scale t^), but 
the level populations relax to their asymptotic value in 
a longer time, two order of magnitude larger than t^,. 
After very large times population inversion arises, with 
Pi > Pq. Indeed, population inversion in the asymptotic 
regime of TLS has been found in 27-2^ under different 
approaches. Within the approximation of (29| it can be 
understood as mediated by virtual transitions to bath os- 
cillator states at energies Eo + nhoj. Here we find that the 
asymptotic regime is difficult to reach, since population 
inversion needs long times (t ^ t^) to emerge when me- 
diated by the bath. We also find that the inversion effect 
and the difference between the finite time and asymptotic 
Pr are enhanced when decreasing temperature (see inset 
of Fig.l) [33. 



(a) Third-level-mediated population inversion: For 
higher fac, the E2 and E3 levels are involved in the dy- 
namics. In this case, the relaxation of the populations Pi 
to the asymptotic regime takes place fast, in a time scale 
similar to t^. In Fig. 2(b) we see that first Pq decreases 
while the populations Pi,P2,Pa increase (through a se- 
ries of LZS transitions, see [llj,|25|), until at a later time 
P2 and P3 decrease transfering population to Pi. An 
asymptotic regime is quickly reached where there is pop- 
ulation inversion. Pi > Po, mediated by E2 and E3, an 
effect that can be observed even in time scales t>tA,. 
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FIG. 3: (color online) Population Pr as a function of the 
resonance frequency of the dc SQUID fip for T = 20mK, 
fdc = 0.50151, fac = 0.00245, uj = 0.003, 7 = 0.001 and 
K = 0.001. Continuous line shows Pr in the stationary state 
and dashed line corresponds to calculations performed for t = 
lOOOr. The blue circles represent the values of A+ and A_. 
The inset shows a detail of Pa as a function of Q.p around the 
n = 9 resonance. See text for details. 

The measurement of the state of the DJFQ is done 
with a read-out dc SQUID, which is inductively coupled 
to the qubit 0, [13, ISH. After considering the measur- 
ing device as a part of the environment, the bath spec- 
tral density can be modeled as Jsh(w) — 7wr2p/[(r2p — 
ijj'^Y + {2nKUjVlpf'], with fip the resonant frequency of 
the SQUID detector, and 2iTKVLp the width of the reso- 
nance at w = Op [23|. Here we calculate the asymptotic 
Pfl solving Eq.(2) but now with an structured bath spec- 
tral density Jsb{(^)- In Fig. 3 we show the results for Pj^ 
as a function of ftp for the same fac of Fig. 2(a). We 
find three regimes in the overall behavior of Pr vs Vlp 
that we describe below. Defining A(/) = Ei{f) — Eo{f), 
A+ = A{fdc + fac) and A_ ^ A{fdc - fac) we get: (1) 
Multiphoton resonances with the hath. For Q,p < A_ the 
population Pr has a strong dependence on Op, displaying 



resonances for Qp « Ticj+ei— eo (with ei, eo the quasiener- 
gies of the Floquet states mostly weighted in the two low- 
est Hq eigenstates). These resonances correspond to the 
maxima at fip of the Jsb{£a — £(3 + nhw) appearing in the 
transition rates in Eq.(2), and have been analyzed in 22 1 
within a perturbative approach. Near these resonances 
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it is possible to fully control the population Pn, going 
from Pji ~ 1 to complete population inversion ~ 0, 
with small changes of flp. In the inset of Fig. 3 we show 
in detail the case near the n = 9 resonance. (2) Bath 
mediated population inversion. In the opposite limit, 
rip > A+, the behavior is similar to the one discussed 
previously in Fig. 2(a), observing population inversion at 
very long times mediated by the bath (note that for large 
flp, the spectral density Jsfc(w) tends to the ohmic one). 
(3) Bath mediated cooling. In the intermediate regime, 
A_ < rip < A+, virtual LZS transitions to an effective 



tions 27H29l| has not been observed in the experiments 



energy level [22, l3l| at Eo + flp < Ei predominate. From 
this effective level it is possible to decay an repopulate 
the ground state, obtaining Pn 1. The resonances at 
rip Ki noj + €i — eo are still observed. A full picture of 
the effect of an structured bath can be observed in Fig. 4, 
which shows a map of the asymptotic population P^ as 
a function of ilp and fac- The above mentioned three 
different regimes for ftp < A_, A_ ^ ^ip < A+ and 
A+ < rip can be distinguished. Also, in Fig. 4 we ob- 
serve the onset of LZS oscillations involving E2 and E3 
at higher fac (corresponding to the second 'diamond' of 
[Tli] ) where the third level mediated population inversion 
discussed in Fig. 2(b) takes place. Here we see that while 
the bath-mediated population inversion mechanism ap- 
pears for rip > A_|_, the third- level-mediated population 
inversion mechanism is almost independent of Op. 




FIG. 4: (color online) Contour map of the stationary pop- 
ulation Pr as a function of the driving amplitude fac and 
the resonance frequency of the dc SQUID Jlp for T — 20mK, 
fdc = 0.50151, uj = 0.003, 7 = 0.001 and k = 0.001. The 
black and green lines represent the positions of A+ and A_. 
Black circle: fac corresponding to Fig. 2(a) and Fig.3; black 
square: fac corresponding to Fig. 2(b). 

In conclusion, by performing a realistic modeling of the 
flux qubit we are able to analyze different scenarios for 
population inversion in strongly driven quantum systems, 
understanding the parameter regimes for their occur- 
rence. One puzzling situation is that the bath-mediated 
population inversion expected to happen in LZS oscilla- 



of [llj. Here we show that this effect is observable only 
when driving the system for very long times, after full 
relaxation with the bath degrees of freedom {t ^ tr). 
The large time scale separation tr S> t</, present in the 
device of 11 ill ( due to the very small gap A = Ei — Eq 
at f = 1/2 [24l. [3^) explains that the predicted popula- 
tion inversion is beyond the experimental time window. 
It will be interesting if experiments could be carried out 
for longer driving times in this device. Population in- 
version in LZS oscillations has been observed in another 
flux qubit device driven at large frequencies [s^, and 
interpreted in [3^ as due to quantum noise, (which is 
consistent with our asymptotic Pr results). Another in- 
teresting point we found is the dependence of the LZS 
oscillations on the frequency rip of the SQUID detec- 
tor. The flux qubits fabricated with Nb junctions as in 
lol 11 1 have typically a small gap and thus they are in 
the regime rip > A+, where bath- mediated population 
inversion takes place. The flux qubits fabricated with Al 
junctions, as in [1, [sH, have a large gap and thus they 
are expected to be in the regime rip < A_. Amplitude 
spectroscopy measurements like ll| but carried out in 
Al flux qubits could give diflterent 'diamond' patterns as 
function of fdc, fac with resonance and LZS intereference 
effects from the oscilator levels at Ei + nhilp. In princi- 
ple, the frequency rip can be varied (in a small range) by 
varying the driving current of the SQUID detector [31j 
or with a variable shunt capacitor. In this case there is 
room to fully control the population Pr and the degree 
of population inversion as a function of rip. 
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